Edge effects in electrostatic calibrations for the measurement of the Casimir force 
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We have performed numerical simulations to evaluate the effect on the capacitance of finite size 
boundaries realistically present in the parallel plane, sphere-plane, and cylinder-plane geometries. 
The potential impact of edge effects in assessing the accuracy of the parameters obtained in the 
electrostatic calibrations of Casimir force experiments is then discussed. 
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1. INTRODUCTION 

The Casimir force [l| has been demonstrated in a va- 
riety of experimental setups and geometries, yet there 
is an ongoing reanalysis of the level of accuracy with 
which it has been determined, which is crucial for as- 
sessing reliable limits to the existence of Yukawa forces 
of gravitational origin predicted by various models 043] ■ 
In performing Casimir force measurements a crucial role 
is played by the related electrostatic force calibrations, 
and mastering all possible systematic effects in the lat- 
ter is mandatory to assess the precision of the former. 
The presence of anomalous exponents in the power-law 
dependence upon distance and the dependence on dis- 
tance of the minimizing potential [j| have been identified 
as a possible source of systematic effects, to be carefully 
scrutinized in each experimental setup . The distance- 
dependence of the minimizing potential has been recently 
modelized theoretically, after a first effort reported in [7| , 
in terms of non-equipotential conducting surfaces due to 
random patterns of patch charges @. Here we report 
results on a further potential source of systematic error 
by studying, by means of numerical simulations using fi- 
nite element analysis, the influence of the unavoidable 
presence of edges on electrostatic calibrations in various 
geometries (for the influence of edge effects on Casimir 
forces in peculiar geometries see 

IS EH)- We limit the 
attention to the three geometries that, apart from the 
crossed-plane investigated in [11], have been extensively 
discussed for measuring Casimir force so far, i.e. the 
parallel plates [H, Ell > the sphere-plane (3-E3] , and the 
recently proposed cylinder-plane |l8( (see Fig. [TJ. We 
discuss the deviations from ideality in all these geome- 
tries through the essential knowledge of the capacitance 
dependence on distance, and its general implications for 
determining the parameters used in the measurement of 
the Casimir force. This is performed under the simplyf- 
ing assumptions that the surfaces are equipotential, i.e. 
by omitting any superimposed effect due to electrostatic 
patches, and by neglecting the tensorial nature of the ca- 
pacitance among all the conducting surfaces realistically 
involved in concrete experimental setups. 



2. PARALLEL PLATE GEOMETRY 



We start our analysis with the parallel plates configu- 
ration since this is the simplest geometry even in terms 
of possible deviations from the ideal, infinite plane case. 
Moreover, besides having an exact expression in the ideal 
case of infinite plates, analytical approximate expressions 
are also available for the capacitance including edge ef- 
fects [3], allowing to obtain reliable numerical bench- 
marks. We consider two identical parallel square plates 
of length L, and boundary conditions set in such a way 
that the two plates are at a constant electric potential dif- 
ference. The total electrostatic potential energy is then 
computed by numerically solving the Laplace equation 
using a dedicated finite element analysis software (COM- 
SOL). We evaluate the total electrostatic energy W e \ by 
summing the electrostatic energy density over a selected 
volume surrounding the two equipotential surfaces. The 
capacitance C can then be calculated from the relation- 
ship C = 2W C \V 2 . The parameters of the mesh (mesh 
size, rate of growth etc.) are carefully chosen and ex- 
tensively tested to ensure the accuracy of the numeri- 
cal results. When the distance between the two plates 
is much smaller than L, the capacitance obtained from 
the numerical simulation agrees within 0.01% with the 
well known analytical formula C pp = e^A/d, where A 
is the surface area of the plates, d their separation, and 
eo the vacuum electric permittivity. However, when the 
distance becomes larger, the value of the capacitance be- 
gins to deviate from the analytical formula. The left plot 
in Fig. [2] shows the capacitance of two parallel square 
plates vs. distance, together with the expected capaci- 
tance for two plates following the infinite surface formula, 
appearing as a dashed line. The data at large distance 
deviate significantly from the dashed line, and the fact 
that they are all above implies that the power-law expo- 
nent is softer, i.e. in between and -1. This deviation 
is ascribed to the finite size of the parallel plates. In 
fact, if only the field lines in the volume delimited by the 
two plates are used for the sum of the electrostatic en- 
ergy, which corresponds to ignoring the contribution of 
the outer region, the capacitance obtained from the nu- 
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FIG. 1: Selected geometries for the numerical study of finite size effects. Parallel plates of finite size (a); a sphere in front of a 
finite size plane (b); a cylinder in front of a plane, with width smaller (c) and larger (d) than the size of the plane. In the latter 
three configurations, a truncated sphere and short and long truncated cylinders, as shown by the dashed lines, have been also 
studied. 



merical simulation would still agree with the analytical 
formula even at the largest explored distances. In order 
to better quantify this deviation we have fitted the ca- 
pacitance curve by progressively removing points at the 
largest distance. As shown in the right plot in Fig. [21 the 
optimal exponent becomes smaller than unity when data 
at large distance are progressively included in the fit. 



3. SPHERE-PLANE GEOMETRY 

In the case of the sphere-plane geometry, the exact 
expression for the capacitance between a sphere and an 
infinite plane is written in terms of a series [20| : 



C sp = AttcqR sinh(a) } — 
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with cosh(a) = 1 + d/R, R is the sphere radius, and 
d the separation distance. In the limit of small sepa- 
rations, d/R << I, an approximate expression for the 
capacitance can be obtained [2fij ]: 
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where < 6 < 1. By neglecting the distance- independent 
terms one obtains an expression often appearing in rela- 
tionshi p to the so-called Proximity Force Approximation 
(PFA) [2l|,|22j in electrostatics. The formulas above both 
assume an infinite plane and a whole sphere. In real ex- 
periments involving microresonators the size of the plane 
is not necessarily large enough to be considered infinite 
and thus edge effects could be present. In some measure- 
ments 17[ the sphere is located close to one end of the 
squared plane, to increase the torque exerted on the un- 
derlying microresonator. Also, in other measurements a 
lens, schematized as a truncated sphere was used in lieu of 
a whole sphere 14,23|. Fig. [3]shows the numerical results 
for configurations taking into account these deviations 
from the idealized well as the curves expected 

from the exact ([T]) and the approximate ^ expressions. 



For the whole sphere, it is worth remarking that neither 
the exact nor the approximate expression can give an ac- 
curate value for the capacitance between a sphere and 
a finite plane, even at small distance. We have checked 
that the idealized case of a sphere and an infinite plane is 
approached by considering square plates of progressively 
larger size. The capacitance in the realistic case still pre- 
serves a logarithmic dependence at small distance, with 
the slopes very close to each other. For the truncated 
sphere, the capacitance is significantly smaller than both 
the exact and the approximate expressions, which is ex- 
pected considering that there is less conducting surface 
available in this case. Moreover, much smaller distance 
is required for the capacitance to be approximated by a 
logarithmic dependence with a slope comparable to the 
one of the whole sphere case. The truncation aspect ra- 
tio has been chosen in close analogy on the case of the 
lenses used in the sphere-plane experiment reported in 
[23j, and it should also be similar to the one used in 
the first modern Casimir force experiment in the sphere- 
plane configuration 14 [. 



4. CYLINDER-PLANE GEOMETRY 

For the cylinder-plane geometry, the expression for the 
capacitance between a cylinder and an infinite plane is 



cp ~ coshT^l + 



d/R) 



(3) 



where L and R are the length and radius of the cylin- 
der, and d the distance. In the limit of small separations, 
d/R « 1, an approximate expression for the capaci- 
tance can be obtained by expanding the inverse hyper- 
bolic cosine function with Puiseux series [24j: 
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Beside the finite size of the plane, the finite length of the 
cylinder may also be a source of deviation from ideality. 
Therefore we have considered two cases, differing in the 
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FIG. 2: (Left) Capacitance vs. distance for two parallel square plates. Dots indicate the results of the numerical simulation, 
while the dashed line is the capacitance expected using the formula C pp = eoA/d. The length of the square plates L was chosen 
to be 8.86 mm. (Right) Optimal exponent e c coming from the best fit of a set of numerical data with C pp = eoA/(d — do) ec vs. 
the distance of the farthest data point used in the fitting. Red dots are obtained by fixing <fo=0, corresponding to an a priori 
knowledge of the absolute distance, while the black squares are obtained if do is considered as a fitting parameter, as usually 
done when analyzing real experimental data with no a priori independent knowledge of the absolute distance between the two 
surfaces. Considering the size of the square plates used in the experiment reported in pl| of 1.1 mm, all separation distances 
in that case are obtained by scaling down the horizontal axis by a factor ~ 65. 



length of the cylinder chosen as smaller or larger than 
the side length of the plane, and the analysis has been 
repeated for truncated cylinders. The capacitance per 
unit length C cp /L e ft vs. distance for different configura- 
tions are shown in Fig. 01 In this figure L e s is defined as 
the minimum between the length of the cylinder and the 
length of the plane along the cylinder axis, i.e. defining 
the overlapping area between the cylinder and the plane. 
As can be seen in Fig. 21 truncated, wider cylinders devi- 
ate less from the exact expression in comparison to whole, 
narrower cylinders. However, longer width means harder 
parallelization and consequently worse achievable mini- 
mum distance. From Fig. 21 a narrow truncated cylinder 
seems to be a good compromise, which is the configura- 
tion used in realistic experiments. 



5. IMPLICATIONS FOR THE ELECTROSTATIC 
CALIBRATIONS IN CASIMIR FORCE 
EXPERIMENTS 

The discussion above shows that when the sizes of the 
samples are finite and the edge effects cannot be ignored, 
the capacitance will not scale with distance through the 
exponent expected from the ideal formula. If the expo- 
nent is incorrectly forced to the value from the ideal for- 
mula in the data fitting procedure, systematic errors will 
be propagated to all the fitting parameters. In fact, the 
optimal exponent is not a constant and depends on dis- 
tance as evident in the plane-plane configuration shown 
in Fig. [21 therefore assuming a constant optimal exponent 
over the entire range of explored distances will generate 
systematic errors in the electrostatic calibrations. 



Let us consider the cylinder-plane geometry as an ex- 
ample. As already mentioned, in a real experiment it is 
hard to determine directly the absolute separation dis- 
tance, and only data fitting allows to obtain the sepa- 
ration if considered as a fitting parameter (see also the 
detailed discussion in [25| on the same definition of dis- 
tance between two macroscopic bodies) . In the ideal case, 
by fitting the capacitance data with a relationship like 

C cp = Co+K c /(d-d Q yc (5) 

with the exponent e c = 0.5 and Co, K c and do as fit- 
ting parameters, the absolute distance can be obtained 
as d a bs = d — do. Here Co represents a constant fit- 
ting parameter due to background, parasitic, distance- 
independent capacitance added to the cylinder-plane ca- 
pacitance, obtainable in the limit of very large gap sepa- 
ration. In the case of our numerical simulation data, we 
can fix Co = 0. If the correct fitting equation is used, the 
value of do obtained should agree within the fitting errors 
with the expected null value. By fitting with Eq. ([5]) the 
simulation data of the narrow truncated cylinder up to 
distance 10 /im, the maximum distance at which electro- 
static calibrations are usually carried out, do = — 18 ± 3 
nm is obtained with the exponent e c is fixed at 0.5, while 
we obtain do = 12 ± 1 nm if e c is left as a free parameter, 
with an optimal value of 0.4849 ± 0.0005. In both cases, 
the value of do from the fitting is different from zero by 
several standard deviations. The smallest distance in the 
data is 400 nm, -18 nm and 12 nm represent errors with 
respect to the smallest gap separation of 4.5 % and 3.0%, 
respectively. Similar analysis can be repeated for the par- 
allel square plates, and the whole and truncated sphere 



4 



• Centered 
■ Lateral 

♦ Truncated 

Exact 

PFA 

IPFA 



c 




1.0x10 

d(m) 



c 

u 

u 



Truncated (w=4 mm) 
Whole (w=4 mm) 
Truncated (w=12 mm) 
Whole (w=12mm) 
Exact 



1.0x10" 

d(m) 



FIG. 3: Capacitance dependence on distance in the sphere- 
plane configurations. The case of a whole sphere located 
above the center of a square plate (black dots) and on the edge 
of the plate (red squares) are shown, as well as the case of a 
truncated sphere (blue diamonds). The radius of the sphere 
is chosen to be 0.15 mm and the length of the square plate 0.5 
mm. The continuous line corresponds to the exact expression 
as in JTJ), and the dashed dot line corresponds to the approx- 
imate expression as in ((2j, named Improved PFA (IPFA), 
while the dashed line is the capacitance expected from the 
PFA method. 



FIG. 4: Simulations in finite-size cylinder-plane configura- 
tions. Capacitance per unit of effective length vs. distance 
for a truncated cylinder (black dots) and a whole cylinder 
(red squares) of width 4 mm in front of a square plane of 
size 10 mm x 28 mm, and similar plots for truncated (blue 
diamonds) and whole (green triangles) cylinders of width 12 
mm. The radius of the cylinder is chosen to be 12 mm. The 
continuous line corresponds to the exact expression as in 
and the approximate expression Q coincides with the exact 
expression in this distance range, confirming that the PFA 
holds to a high degree of accuracy up to the largest explored 
distances. 



configuration by fitting the data respectively with 

c pp = c + K c /(d-d y% 

C sp = C Q + K c \n{d-d Q ). 



(6) 



The results are shown in Table I, where the values of 
do are all significantly different from the expected null 
value. To show that this is not a fitting artifact, the same 
fitting procedures have been repeated with simulation 
data at much smaller distance where the edge effects can 
be safely ignored. In this check the values of do agree 
with the expected value of zero within the associated 
error from the fitting procedure. Therefore it is evident 
that the presence of edge effects may result in significant 
systematic errors in the best fit, affecting the accuracy of 
crucial parameters such as the absolute distance. 

With regard to Casimir experiments, other crucial pa- 
rameters of the experiment beside the absolute distance 
between the two surfaces arc determined through the 
electrostatic calibrations, for instance the effective mass 
m e s of the resonator or, equivalently, its stiffness. The 
calibrations are usually done at relatively large distance 
to ensure that the Casimir force can be neglected. How- 
ever, in this regime the edge effects become more signifi- 
cant leading to systematic errors in the fitting procedures 
as shown above. However, in most Casimir experiments 
what is measured in the calibrations, rather than the 
capacitance, is the electrostatic force between the two 
surfaces or, if a resonator is used for dynamical measure- 
ments, its square frequency shift [26j . 

The relationship between these two obscrvables and 



the capacitance is in general given by 
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Then the absolute value of the capacitance in itself is ir- 
relevant for force or square frequency shift measurements 
and, due to the first and second derivatives, the distance 
dependence of the force or frequency shift signal may end 
up being less sensitive to the edge effects. 

To check the influence of the edge effects on force mea- 
surements, numerical differentiation has been carried out, 
using the Lagrange three-point formula [13] , on the sim- 
ulation capacitance data. The numerical differentiation 
method works better when the data points are equally 
spaced, therefore in the case of the sphere-plane con- 
figuration the numerical differentiation is performed in 
the semilogarithmic domain, while for plane-plane and 
cylinder-plane configurations in the logarithmic domain. 
Then we fit the force vs. distance data with the corre- 
sponding equations for each geometry, obtaining the re- 
sults shown in Table II. It seems that in terms of the op- 
timal exponent, the edge effects are less significant when 
the force data are considered, apart from the case of the 
truncated sphere. More specifically, the results for the 
truncated sphere can be understood as if it behaves like 
a sphere only at small distance, instead resembles a plane 
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Parallel plates 


Sphere I 


Sphere II 


Cylinder 


do(nm) 

do/ dmin 


-90 ± 20 
1.8% 


13 ±3 
6.5% 


144 ± 13 
72% 


-18 ±3 
4.5% 


d' (nm) 

dQ / drain 


140 ± 20 
2.8% 
0.977 ±0.002 


N/A 
N/A 
N/A 


N/A 
N/A 
N/A 


12 ± 1 
3.0% 
0.4849 ± 0.0005 



TABLE I: Parameters of the best fit obtained for the capac- 
itance data in the case of the four configurations considered 
in the text with Sphere I being whole sphere and Sphere II 
truncated sphere. For parallel plates and cylinder-plane con- 
figurations, the top two rows are the parameters obtained 
when the exponent e c is fixed at 0.5 and 1, and the bottom 
three rows when e c is left as a free parameter. For the sphere- 
plane configuration, due to the logarithmic dependence of the 
capacitance scales with distance, no exponent can be intro- 
duced. 



when considered at large distance from the plane. This 
result is at variance with respect to the behavior actu- 
ally observed in since a square frequency shift power 
law dependence on distance with an exponent less than 2 
(therefore corresponding to a force exponent ef < 1) has 
been observed uniformly over the entire range of explored 
distances. Therefore the origin of the anomaly observed 
in [H[ cannot be ascribed to edge effects. Analogous con- 
siderations can be repeated for the first derivative of 
the force whenever square frequency shifts are measured. 
The influence of the edge effects on the values of fitting 
parameters do and K{ in force measurement is also less 
significant at relatively small distance in the case of par- 
allel plates, truncated cylinder and whole sphere, while in 
the case of the truncated sphere the error is much more 
significant. However, although the edge effects are not 
as strong in force measurement as in capacitance mea- 
surement, related systematic errors are still present, es- 
pecially at large distance. Therefore the edge effects need 
to be considered and carefully examined if high precision 
measurements are required, and the robustness of the fit- 
ting procedure with respect to the choice of the distance 
range becomes crucial for the parameters, in particular 
Kf and do- 



6. CONCLUSION 

We have analyzed the effect of finite size in the pre- 
dictions for the electrostatic calibrations of geometries of 
interest for measuring the Casimir force. One relevant 
consideration arising, or better quantitatively confirmed, 
from this numerical study is that there is a trade-off in 
the choice of the distance range to perform electrostatic 
calibrations. If this range is including too close distances, 
roughness corrections and the presence of the attractive 
Casimir force may result in systematic effects, while if 
large distance data are also considered one must include 



edge effects in the theoretical modclization of the exper- 
imental setup for a targeted level of accuracy of the elec- 
trostatic calibrations. If geometrical configurations man- 
ifestly differing from the idealized cases are also consid- 
ered, such as a small portion of a sphere or a truncated 
cylinder, the deviations in the capacitance are even more 
pronounced and a full numerical evaluation of the elec- 
trostatics is necessary even at relatively moderate levels 
of desired accuracy. We have also shown that force (gra- 
dient force) measurements, based on the first (second) 
derivative with respect to distance of the capacitance, are 
in general more robust than absolute capacitance mea- 
surement in itself with respect to the geometrical devia- 
tions from ideality, confirm ing from another perspective a 
point recently discussed in [28| . Finally, we want to point 
out that detecting edge effects provides also a simple and 
reliable method to assess the precision of the electrostatic 
calibrations, for instance determining the minimum sep- 
aration gap above which deviations attributable to finite 
edges are effectively observed. 
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Af(Nm 'V ) 
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